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^ Abstract 

o 

In this paper we introduce a completely continuous and time-variate model of the evolu- 
tion of market limit orders based on the existence, uniqueness, and regularity of the solutions 
to a type of stochastic partial differential equations obtained in [6 . In contrary to several 
models proposed and researched in literature, this model provides complete continuity in 
both time and price inherited from the stochastic PDE, and thus is particularly suitable for 
Ch the cases where transactions happen in an extremely fast pace, such as those delivered by 

1^ high frequency traders (HFT's). 

I ^ i ' We first elaborate the precise definition of the model with its associated parameters, 

^ and show its existence and uniqueness from the related mathematical results given a fixed 

^ set of parameters. Then we statistically derive parameter estimation schemes of the model 

using maximum likelihood and least mean-square-errors estimation methods under certain 
criteria such as AIC to accommodate to variant number of parameters . Finally as a typical 
economics and finance use case of the model we settle the investment optimization problem 
in both static and dynamic sense by analysing the stochastic (Ito) evolution of the utility 
OA function of an investor or trader who takes the model and its parameters as exogenous. Two 

• • theorems are proved which provide criteria for determining the best (limit) price and time 

point to make the transaction. 
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1 Introduction 

In literature there have been researches on the modeling of market limit orders and their execu- 
tion such as Uj and [5], most of which are based on discrete settings, for instance, models based 
on Poisson processes and/or queuing theory. However, because of rapid technological evolution 
which brings about ultra- fast microprocessors and hardware, trading behaviors and patterns in- 
volved with a large amount of limit order creation, transaction, and cancellation within a short 
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Figure 1: An Illustration of Moving Boundary PDE's 

period of time, such as high frequency trading (HFT), have become quite popular and tend to 
have a heavy impact on the mechanisms of price discovery and formulation. Consequently, a 
continuous and dynamic model of the evolution of limit orders in a particular market and their 
dynamics is proposed, which is described by a type of stochastic partial differential equations, 
Stefan equations, with a number of model parameters to be estimated given a real dataset. 

Such type of equations of u{t, x) describes the behavior of a system that consists of two 
phases, as illustrated in Figure [l| where f3{t) is a moving boundary which is part of the solution 
and must be solved simultaneously with u{t, x). As can be seen in Figure [l| in the region to the 
left of the moving boundary (namely, the set {{t, x) £ [0, oo) x M : x < /3(t)}) u is constantly set 
to 0; on the right side of the boundary (the set {{t, x) £ [0, oo) x M : a; > /3(t)}) u is described by 
a PDE of the general form Lu + 6 = where L is a predefined second-order differential operator. 

For a moving boundary PDE problem, in addition to the regular boundary condition such as 
the Dirichlet condition u{t, (i{t)) = 0, there is always an extra boundary condition that describes 
the dynamics at the moving boundary, for instance, the Stefan boundary condition 

rill 

-{t,m+) = pm- (1) 

The type of moving boundary PDE's we shall use in modeling of market limit order evo- 
lution is the Stefan problems, where L is a heat or parabolic operator (for instance, L := 
—d/dt + d"^ /dx"^) with the Stefan boundary condition. Such type of problems has a variety of 
applications. For instance, in physics, they model the phenomena such as ice melting with the 
Stefan condition describing the heat balance at the interface (the moving boundary, see [2]); in 
finance they model the valuation of American options with the PDE derived from the Black- 
Scholes formula and the moving boundary describing the early exercise price boundary (see 
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Lemma 7.8, Chapter 2, [3]). As it turns out in this paper such equations are also perfectly suit- 
able for modeling the dynamics of high frequency limit order transactions and their evolution, 
providing statistical parameter estimation methods, and formulating and solving the investment 
optimization problem based on a typical utility function. Note that [6j provides the mathe- 
matical foundation and develops essential techniques that enable us to justify the existence and 
uniqueness of such type of equations, while obtaining its essential regularities. 

In Section [2] we model the evolution of limit orders in a particular market by Stefan equations 
according to the following facts, and show the existence of such model based on the mathematical 



(1) Limit orders are placed, cancelled, and executed in a manner where jitters tend to be 
rapidly smoothed out, which is why we have a Laplacian; 

(2) The change of the mid-price is driven by the intensity of interaction between ask and bid 
orders around the mid-price; 

(3) The randomness comes from the constant creation, cancellation, and execution of limit 
orders; its intensity varies at different (limit) prices, and tends to vanish as the price goes 
far beyond the mid. 

In Section |3] we study the methods to estimate the model parameters based on a given limit 
order dataset and derive the statistics such as a Maximum-Likelihood Estimator (MLE) and 
an estimator that minimizes the Mean-Square Errors (MSE), both under AIC (see [Tj) since 
the number of parameters (or dimension) is to be estimated itself. Under certain simplified (or 
degenerated) circumstances an explicit expression of an MLE estimator is also derived. 

In Section|4]we study the investment optimization problem derived from the model, formulate 
and analyze the static and dynamic properties of the utility function of an investor who takes the 
model and its parameters as exogenous, which finally enable us to find the optimal limit-price-to- 
buy (or equivalently, amount-to-buy) of the asset and the corresponding amount of consumption, 
given a fixed amount of wealth and at a given time. Both static and dynamic analysis are studied 
to obtain the optimality of the investment via the limit order model, and two theorems are given 
respectively for those two types of analysis as the criteria to test for optimality from the model 
dynamics using Ito's Formula, which also have intuitive interpretations. 



We model the evolution of the limit orders of a particular asset in the market. Suppose we 
work within the time interval [0,T], and S denotes the natural log of the price so 5 G M. At a 



results obtained in 
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S*{T) 




Figure 2: The evolution model of the limit order book 

particular time point t £ [0,T], suppose the volume of the ask (resp. bid) limit orders on the 
limit order book from S to S + dS is VAit,S)dS (resp. Vsit, S)dS), and suppose the natural 
log of the mid-price is S*{t). Fix a probability space (0, 3~, P), We then have the following fairly 
natural but important model assumptions: 

(1) all matched orders are executed immediately since all major trading centers today have 
been computerized, then we have for all {S,t) such that S > S*{t), VB{S,t) = 0; similarly 
for ah {S,t) such that S < S*{t), VAiS,t) = 0; 

(2) since jitters have a tendency to be rapidly smoothed out, dVA/dt contains a component 
CKA'^Kl where oa is a positive constant and A := d'^/dS'^ is the Laplacian; the same is 
true for Vb^S^I) with as', 

(3) the change of mid-price is driven by the "strength" of the ask and bid orders placed around 
the mid-price, which implies we have a Stefan- type condition 



dS% , 



(2) 



where p is a constant; 



(4) for S > S* (t) the ask orders are placed in a stochastic manner, hence another component of 
dVA/dt is <yA{\S—S* {t)\)j^^ where a a is a regular scaling function (defined in Chapter ??) 
and W :Vt X [0, T] x M — )• M is a standard 2-dimensional Brownian sheet; the same is true 
for Vb with S < S*{t) and as- 
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The model is illustrated in Figure [2j 

In sum, for a set of given model parameters uqa, uqb, o:a, oBt^ta, o'b, P, the model has 



OVa 



+ ctaHs - s*m—,ys > s*{t), 



dt dS^ 
VA{t,S) = 0,V5 < S*{t), 

Va{0,S) = uoa{S), 

d'VB 



dVB 



OiB- 



dt ^'^ 
Vsit^S) = 0,V5 > 5*(t), 

dVA 



+ aB{\S-S*{t)\) — ,yS<S*{t) 



dtdS' 



(3) 



dS\ , 



dS 



{t.S*{t)+) + ^-^{t,S*{t)- 



Theorem 2.1 The solution Va,Vb, S* to model ([sj) exists and is unique for <t < t where r 
is a well-defined stopping time. 

Proof We roughly follow the same procedure to show the existence and uniqueness as in [Hj. 
First we make the following transformation: VA{t,S) := VA{t, S* (t) + S),uoAiS) := uoa{S* (0) + 
S),VB{t, S) := VB^t, S*{t) — S),uqb{S) := no_B(5'*(0) — S). Then the same argument as in 
Theorem 3.2, [6j shows that there exists unique (SajPb that is the limit of the iteration 

pPAit) = J ^it,0,y)uoA{y)dy + J J ^{t - s,0,y)aA{y)Wpn{dyds), 

pl3Ut) = J^ ^{t,0,y)uoB{y)dy + ^{t - s,0,y)aBiy)Wp^idyds), 

p^{t) = pm-m)- 

Define S*{t) := lim„/3"(t) which is also unique. Then using the same argument as in Theorem 
3.2, [6j, we have that and Vb exist and are unique. Also, using Kolmogorov's Continuity 
Theorem as in Lemma 3.1.1, we have that almost surely, 

/^/3A(t) = ^(t,5*(t)+), 
P/3B(t) = -^(t,5*(t)-), 

where < t < r := limi^ooT^, and := inf{t > : \pA{t)\ > L or \$B(.t)\ > L}. Therefore 
the Stefan boundary condition ([2]) holds. □ 
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3 Parameter Estimation 



In this section a statistical method based on AIC is developed to estimate the parameters of 
model ^ given a real dataset of the limit order book and the mid-price of a particular asset 
within a certain period of time. 

Suppose the dataset consists of 3 matrices, two T x N matrices D^, Db for the volumes of 
the ask and bid orders, and a T x 1 matrix P for the mid-price. The sampling steps for time 
and price are At and Ajy, that is, the dataset is from time to Aj^T, and the t-th row of Da 
stores the volumes of the ask orders from limit price P[t\ to P[t\ + A^rA^, while the t-ih. row of 
Db stores the volumes of the bid orders from limit price P[t\ — AjvA to P\t\, where [•] denotes 
the vector subscription. 

Our goal in this section is to develop an algorithm to numerically compute or even find an 
explicit expression of the statistics served as the appropriate estimators of the model parameters 
under maximum-likelihood, minimum mean-square errors, and AIC. The method is decomposed 
into two parts, first parameters a^i, a^, cta, are estimated using maximum-likelihood ap- 
proach combined with AIC based on the limit order book data Da and Db, and then uqa, uob, P 
are estimated using least mean-square-errors combined with AIC based on the whole dataset 
Da,Db,P. 

3.1 AIC/MLE Estimation of a^, a^, a^, as 

To estimate a a and cjb, we further assume that 



where for i = A, B , pi is a, polynomial with deg pi > 1, so that cjj is a regular scaling function, 
and lim^^-i^oo o'i(x) = 0, because there tends to be very little randomness as the limit price goes 
far beyond the mid- price. Then the goal is to estimate aA,OiB and the degrees of pa,Pb and 
their coefficients. Since the number of parameters 4 -|- degpA + degp^ is also to be estimated, 
an AlC-based approach is used. 

Suppose i = A, B, and Pi{x) = Yl'j=oPij^'' '^ith di > 1. To estimate ai,pio, . . . ,Pidi, we use a 



finite-difference Euler approximation scheme similar in Section 4, [6j. For each integer n, define 
the difference operator V : M" — )■ M"~^ by 



Denote (Va) [i] by Vcj. For a matrix with row label R and column label C denote as 
the difference operator on column vectors and Vc as that on row vectors. Then for each 



V ((ai, . . . , a„)^) = (a2 - ai, . . . , 
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i = A,B,t = l,...,T -l,S = l,...,N -2, 



+ ai{SA 



At An ' "'^"-^^AtAjv 

where 

Ct,s := W {[ATt, Arit + 1)] x [S*{ATt) + A^S, S*{ATt) + An{S + 1)]) 

are independent identically distributed Gaussian random variables with mean and variance 
AtAat. Fix di = deg Pi, then the log likelihood function £i{ai,pio, . . . ,Pidi) satisfies 



6,5 



T-1 N-2 



-2ii{ai,pio, . . . ,Pidi) = J2^ 



t=l 5=1 



VtDi[t,S]-ai 



VlDi[t,S] 



(4) 

When di is fixed, the first order condition is a cubic formula of the variables, which shall be 
solved using numerical methods. If the model can be degenerated so that = aoi is a known 
constant, then we can indeed solve for optimal Pio, ■ ■ ■ ,Pidi- 

Pi = -A,"^bi, 

where Pj = {pio, . . . ,Pidi)'^ , Aj is a (d, + l)-square matrix with its {m,n) element being 



T-l N-2 

EE 

i=l 5=1 



VtDi[t,S]-ao^[ ^]vlDi[t,S\ 



A 



N 



m+71— 3.2 ^m+n 



N 



and bj is a ((ij + 1) X 1 matrix with its m-th element being 

T-l Af-2 



EE 

t=l 5=1 



VtDi[t,S] - Qoi 



A2 



S 



2m— 3.2 A 2m 



Therefore, under AIC (see [Ij), the set of optimal estimators ai,di,pio, . . . ,pidi minimizes 



AICi := 2di - 2ei{ai,pio, . . . ,PidJ- 
The algorithm to minimize AlCj: 

(1) m := — oo, di := 0; 

(2) For di G {1,...,10}: 

(a) Numerically solve the first order condition of Q for optimal ai,piQ, . . . ,Pid^; 

(b) Calculate a :=AIC(dj) based on ([5]); 

(c) If m > a then m := a, di := di; 



(5) 



(3) The final di is optimal within {1, . . . , 10} with AIC m. 



3.2 AIC/Leeist-MSE Estimation of uqa,uqb,P 

We continue to assume i = A,B. Similarly we need further assumptions about the initial 
condition uoi. Assume 

Uoi{x) = xqi{x) exp{--fix) 

where qi is a polynomial and 7^ > is a parameter, so that uoi has exponential decay at infinity. 
Then Cj := deg qi and its coefficients qiQ, . . . ,qic- arc also a model parameter that needs to be 
determined under AIC. Now, if we take the expectation conditional on S* on both sides of the 
evolution equation of Vi, we get 

roo rt roo 

E[Vi{t,S)\S*]= p.it,S,y)uoi{y)dy+ S*{s) q{t - s, S,y)E[Vi{s,y)\S*]dyds. 
Jo Jo Jo 

This implies that E[l^|S'*] satisfies a deterministic Stefan-type PDE, which means that although 

we are generally unable to find an explicit expression of the solution, we can find a sufficiently 

accurate numerical solution. For a given dataset D^, Db, P, the goal in this part is thus to 

develop an algorithm to minimize the mean-square errors against E[Vi|<S'*] and the mid-price 

function (moving boundary) determined by it under AIC. 

Given a set of parameters ca,cb, qAO, ■ ■ ■ , QAca > QBO, ■ ■ ■ 1 QBcb ilA-ilB-, P-, denote the solution 

to the deterministic Stefan- type PDE (i.e., when W = Q)hy Va and Vb, which can be obtained 

numerically. Then the variances (errors) come from two parts: those from the evolution of the 

limit order book, and those from the evolution of the mid-price. Let ^0 be a predefined constant 

which is the weight of importance of how well mid-price fits against how well the limit order 

book fits. Namely, we shall minimize the weighted MSE 



where 



and 



MSE := MSEi + ^oMSEa 



T N 
ie{A,B} t=i s=i 



MSE2 -.^l^Yl 



^ 'pP(t) yA{ATt,s)^VB{ATt,sy '' 



T 

t=l 

Finally, similar to the AIC approach in the previous part, the goal is to minimize 

CA + CB + MSE, 

which can be done by traversing all possible values of the model parameters, numerically find- 
ing Va and Vb, computing the corresponding MSE's, and finally finding the recorded optimal 
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Figure 3: Investment Optimization via Limit Orders 

parameters. Note that unlike in previous part, because of the nonhnearity of Va and Vb, it is 
difficult to find a direct way to (numerically) solve for optimal coefficients for qi, as opposed to 
Pi- 

4 Investment Optimization 

Using the model combined with a set of optimal parameters adjusted for a given dataset, an 
investor can optimize his or her allocation of the amount of investment in the asset against 
consumption within a given amount of wealth. To analyze the investment behavior under our 
model, we assume the investor can only buy a single asset by making limit order transactions, 
that is, fulfilling ask orders placed by other sellers. This scenario is typical when an investor 
is avert to the volatility of price changes in market order transactions, or there lacks sufficient 
liquidity of the asset but the investor still has a strong motivation to make purchases. 

Let U : Mi]_ — )• M be the utility function of the investor, which is monotone increasing and 
concave down on the second and third parameters. Specifically, U{t,Lt,Ct) denotes the amount 
of happiness the investor obtains at time t if having the amount of asset Lt and consumption 
Ct- Since the investor would choose to pay the lowest possible price to buy an amount of asset, 
the amount of asset Lt and the total cost to buy such amount is completely determined by the 
highest limit price S*{t) + B the investor is willing to pay, which is illustrated in Figure [sj 

Taking the model (and its parameters) as exogenous, and assuming the time t wealth is 
Wt, we then have the following optimization problem at time t (where the subscription t in the 
expectation denotes it is conditional to all information up to time t): 

maxEt [U {t, Lt{B),Ct)] 
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where 

rS*{t}+B 



Lt{B)= / VA{t,S)dS 

Js*H) 



/■S*{t)+B 

Wt = Ct+ / SVA{t,S)dS. 
Js*(t) 



IS*(t) 

subject to the budget constraint 

rS*{t)+B 
is*{t) 

When t is fixed, the optimization of U with respect to the choice of B is done through a 
static analysis; if we consider the full model, that is, both t and B vary, we come up with a 
dynamic analysis. 

4.1 Static Analysis 

Let time t be fixed. Denote the partial derivatives of U by Ut,UL,Uc- Then we compute 

^{t, Lt{B),Ct) = UUt, Lt{B),Ct)VA{t, S*it) + B) 

- Uc{t, Lt{B), Ct)iS*{t) + B)VA{t, S*it) + B). 
Therefore the optimal B = B* satisfies the first order condition 

UL{t,Lt{B*),Ct) 



This means 



Theorem 4.1 In static optimization where the time t is fixed, the optimal highest limit price to 
buy the asset is equal to the ratio of the marginal utility of amount of asset to that of consumption. 

Proof As above. □ 

Note that S* (t) + B* can be seen as the highest amount of money the investor is willing to pay 
to substitute 1 unit of asset with the same amount of consumption. 

4.2 Dynamic Analysis 

Now consider that t also varies, and we are interested in Ei[di7] from t to t + dt. Indeed, if we 
have Et[dU{t,Lt(B*),Ct)] > 0, then even if B = B* reaches its static optimality, the investor 
would still wait for an amount of time to maximize the utility. 
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First we consider the evolution of Lt{B). Fixing B and plugging in cIVa, we have formally 



dLt{B) = a A 



^{t,S*{t) + B)--^{t,S*{t)) 



dS 

S*{t)+B 



+ 



t /■5*(t)+B 



dS 

a{S - S*{t))W{dSdt) 



dt 



a{S - S*it))W{dSdt) 



JS*{t) 



S*{t)dt. 



Similarly, 

dCt = OA 
+ 

+ 



VA{t,S*{t)+B) - iS*it) + B)^it,S*{t) + B) + S*it)^{t,S*it)) 

S*it)+B 

Sa{S - S*{t))W{dSdt) 

5*(i) 

t i-S*it)+B 

/ Sa{S - S*{t))W{dSdt) S*{t)dt. 

Js'-{t) 



dt 



Then by Ito's formula, and noting that S*{t) is part of the information at t, we get 



Et[dU] = Utdt + UiaA 
+ UcaA 



^^^it,S*it) + B)-^it,S*it)) 



dS 



dS 



dt + -ULL / a^{S)dSdt 
^ Jo 



dt 



+ \ucc r {S*{t) + Sfa^{S)dSdt. 
^ Jo 



Since for a fixed t we are only interested at B* in Theorem 4.1 and S* + B* = Ul/Uc, we have 



dU 
~dt 



it,Lt{B*),Ct 



Ut + UcaA 



VA{t,S*{t) + B*) - B*^{t,S*{t)) 



+ \ull r cJ^{S)dS+\ucc r {S*{t) + Sfa''{S)dS. 
^ Jo ^ Jo 



(6) 



If we denote the absolute risk aversions of the investor with respect to the amount of asset 
and the consumption by 

Ull 



rc :-- 



Ul ' 
Ucc 
Uc ' 
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then ([g]) is equivalent to 

^{t,Lt{B*),Ct) =Ut + UcS^aA VA{t, S* (t) + B*) - B*^{t, S* (t)) 

-IvLiS* + B*) Ij a\S)dS - \tc IJ {S* + Sfa\S)dS^ . 

(7) 

The investor would then use the dataset Da to compute [^] and see the expected change 
of U from t io t + dt dX B* . Specifically, we have the following observations from Equation Q 

0: 



or 



(1) the larger the quantity 



is, the more possible it is for U to increase; 

(2) the greater absolute risk aversions and rc the investor has, the more possible it is for 
U to decrease; 

(3) the larger B* is, the more risk U is exposed to for a decrease, since both/g'' a^{S)dSimd 

{S* + S)'^a'^{S)dS are increasing functions of B* . 

Also, if we make two natural assumptions about U, 

(a) the utility function is discounted in time, that is, Ut < 0, 

(b) the investor is risk avert or risk neutral, that is, Ull ^ 0, Ucc ^ 0, 
then we have 



Theorem 4.2 Suppose B* is the statically optimal choice at time t as in Theorem 4-1 ■ If 

.OVa, 



VAit,S*{t)+B*)<B*-^it,S*{t)), 



(8) 



then 



Proof As above. □ 



E. 



dU 

~dt 



{t,Lt{B*),Ct) 



< 0. 



Note that ([s]) in Theorem 4.2 has an intuitive illustration in Figure |4j The quantity 
VA{t,S*{t) + B*)/B* is the slope of the colored lines between the mid-price point {S*{t),0) 
and the optimal point {S* + B* , VA{t, S* + B*)). 
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Figure 4: An Illustration of Theorem 4.2 



* If it is less than the boundary derivative dVA{t, S*{t))/dS as shown in curve Vai, then the 
utility at next instant t + dt is expected to drop, as the volume of ask orders at lower limit 
prices between S* and S* + B* tends to fall, so it is best to make the purchase at t rather 
than wait. 

* On the other hand, as shown in curve Va2, the volume at lower prices tends to rise, and 
it might be wise to wait for a future time to make the purchase. In this case the investor 
needs to evaluate other factors such as time discount Ut and risk premium terms as well. 



5 Summary 

In this chapter we summarize the main results obtained in the previous sections. Throughout 
this chapter we fix a probability space {Q, 3", P), and suppose TyiOxRxM— )-Misa standard 
2-dimensional Brownian sheet. 
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5.1 Modeling 

For a set of given model parameters uqa, uoB,OiA,0!B,o'A,(^Bi P, the model 



dVA 

at 



dt 



exists and is unique. 



VAit,s) = o,ys<s*{t), 

Va{0,S) = uoa{S), 



O.B- 



dt ""^ ds'^ 

VB{t,S) = Q,^S>S*{t), 
Vb(0,5) = uob(5), 

dVA 



+ aB(\S-S*m — ,^S<S*{t), 



'dtds' 



g^ {t,S*{t)+) + ^{t,S*{t)-) 



(9) 



5.2 Pctrameter Estimation 

Suppose P, Da, Db is a given limit order dataset. For i = A,B, suppose 

^1.6 



ai(x) :-- 



,i = A,B 



l + xpi{x)'' 

where pi = Yl'j=oPij^'' ^ polynomial with degpj > 1. Then the first step is to minimize for 
i = A,B, 

AlCj := 2di - 2£i{ai,Pio, . . . ,Pidi) 



T-l N-2 

= EE 

t=l 5=1 



VtDi[t,S]-ai[^)vlDi[t,S] 



(10) 



Suppose also 

uoiix) = xqi{x) cxp{--fix) 

where qi = Yl'j=o 9^^"' ^ polynomial and 7j > is a parameter. Then the second step is for a 
predefined weight to minimize 



CA + CB + MSEi + ^oMSEs 



where 



T N 



MSEi 



2TN 



ie{A,B} t=l 5=1 



MSE. 



2 ■= 



1 ^ 
t=l 



pP(t) 



VA{ATt,S) + VB{ATt,S) 

At 
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and Vi{t,x) is the numerical solution by setting W = 0. 



5.3 Optimization 

Suppose B* (t) is the statically optimal price at time t with respect to the utility function U and 
wealth Wt. Then 



S*{t) + B*{t) 



UL{t,Lt{B*{t)),Ct) 



Uc{t,Lt{B*{t)),Cty 

In other words, the optimal highest limit price to buy the asset is equal to the ratio of the 
marginal utility of amount of asset to that of consumption. 

Also, let rL, rc be the absolute risk aversions with respect to L, C, then we have 



dU 
~dt 



{t,Lt{B*),Ct) 



Ut + Uc< a A 



VAit,S*it) + B*) - B*^it,S*{t)) 
-^rL{S* + B*)£ a\S)dS-^rc£ (S* + Sfa\S)ds\ . 



In particular, if 



then 



B* 

dU 
Itt 



as 



{t,Lt{B*),Ct) 



(11) 
(12) 



< 0. 
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